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(S) Verfahren und Anordnung zum Entwurf eines technischen Systems 

® Zur Losung der Aufgabe wird ein Verfahren zum Ent- 
wurf eines technischen Systems angegeben, bei dem 
MeBdaten eines vorgegebenen Systems anhand eines Er- 
satzmodells beschrieben werden. Es wird eine MaSzahl 
fur die Gute des Ersatzmodells ermittelt, indem die Mefc- 
daten des vorgegebenen Systems mit durch das Ersatz- 
modell bestimmten Daten verglichen werden. Anhand 
der Ma&zahl fur die Gute wird das Ersatzmodell dahinge- 
hend angepafct, daS es eine moglichst hohe Gute auf- 
weist. Das hinsichtlich seiner Gute angepafcte Ersatzmo- 
dell wird zum Entwurf des technischen Systems einge- 
setzt. 
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Beschreibung 

Die Erfindung betrifft ein Verfahren und eine Anordnung zum Entwurf eines technischen Systems. 

Ein Systemverhalten eines technischen Systems, z. B. einer verfahrenstechnischen Anlage oder eines Systems der 
5 GroBindustrie, hangt von zahlreichen Parametern ab. Im Rahmen eines Entwurfs solche eines Systems, also insbeson- 
dere beim Neuentwurf oder bei der Anpassung bzw. Einstellung eines bereits bestehenden Systems, sind Vbrbedingun- 
gen, z. B. hinsichtlich der Wirtschaftlichkeit oder der Umweltbelastung des Systems, einzuhalten. Jede Vbrbedingung 
wird als eine Zielfunktion formuliert, die zu optimieren im Hinblick auf die anderen Zielfunktionen allgemeines Bestre- 
ben ist. 

to Die Aufgabe der Erfindung besteht darin, den Entwurf eines technischen Systems anhand von MeBdaten eines vorge- 
gebenen Systems zu ermoglichen. Gerade im Hinblick auf eine Optimierung des bestehenden Systems bzw auf einen op- 
timierten Neuentwurf eines Systems ist eine derartige Nutzung bekannter MeBdaten von groBer Bedeutung. 

Diese Aufgabe wird gemaB den Merkmalen der unabhangigen Patentanspriiche gelost. Weiterbildungen der Erfindung 
ergeben sich auch aus den abhangigen Anspriichen. 

15 Zur Losung der Aufgabe wird ein Verfahren zum Entwurf eines technischen Systems angegeben, bei dem MeBdaten 
eines vorgegebenen Systems anhand eines Ersatzmodells beschrieben werden. Es wird eine MaBzahl fur die Giite des Er- 
satzmodells ermittelt, indem die MeBdaten des vorgegebenen Systems mit durch das Ersatzmodell bestimmten Daten 
verglichen werden. Anhand der MaBzahl fur die Giite wird das Ersatzmodell dahingehend angepaBt, daB es eine mog- 
lichst hohe Giite aufweist. Das hinsichtlich seiner Giite angepaBte Ersatzmodell wird zum Entwurf des technischen Sy- 

20 stems eingesetzt. 

Die aus vielen unterschiedlichen realisierten Systemen vorliegenden MeBdaten werden zur Beschreibung des Ersatz- 
modells eingesetzt. Mit dem Ersatzmodell wird versucht, das vorgegebene System moglichst gut nachzubilden. Die 
MaBzahl fur die Giite der Nachbildung wird ermittelt, indem die realen MeBdaten mit den Daten, die anhand des Ersatz- 
modells gewonnen werden, verglichen werden. Eine groBe Differenz zwischen den MeBdaten und den Daten des Ersatz- 

25 modells entspricht einer schlechten Giite, also einer schlechten Abbildung des vorgegebenen Systems in das Ersatzmo- 
dell. Durch die MaBzahl fur die Giite wird das Ersatzmodell dahingehend angepaBt, daB die Giite selbst moglichst hoch 
wird und somit das Ersatzmodell das vorgegebene System moglichst gut beschreibt. Das so gewonnene Ersatzmodell mit 
hoher Giite wird zum Entwurf des technischen Systems eingesetzt. 

Unter Entwurf werden allgemein verstanden sowohl der Neuentwurf eines technischen Systems als auch die Anpas- 

30 sung bzw. Adaption eines bereits vorhandenen technischen Systems. Eine Weiterbildung besteht darin, daB das Ersatz- 
modell ein Regressionsmodell ist. Das Regressionsmodell geht von der Beschreibung 

Vi = fp(xO + Gi 

35 aus, wobei 

(yi,Xj) vorgegebene Wertepaare, 

fp eine Funktion, die von einem Parameter P abhangt und 

ei einen Fehler 

bezeichnen. 

40 Nun ist der Fehler (als Funktion von J3) zu minimieren: 

£e? =<p( P ). 
i = l 

45 

Geht man von foigendem Beispiel 
Y = p 0 + PiX + P2X 2 + e 

50 aus, so ist der funktionale Zusammenhang quadratischer Ordnung, das Regressionsmodell (Funktion, abhangig von P) 
hingegen ist linear. 

Die Giite kann in einer anderen Weiterbildung anhand einer quadratischen Abweichung der MeBdaten von den durch 
das Ersatzmodell bestimmten Daten ermittelt werden. Die Anpassung des Ersatzmodells erfolgt durch Minimierung der 
quadratischen Abweichung. 

55 Eine Ausgestaltung besteht darin, daB die MeBdaten nach ihrer Giite, bezogen auf deren Abweichung von den durch 
das Ersatzmodeii bestimmten Daten, sortiert werden und eine vorgegebene Anzahl von n% schlechtesten MeBdaten aus- 
sortiert werden. Es wird also eine Giite fur jedes MeBdatum bestimmt, wobei die Menge der MeBdaten, vorzugsweise in 
Form einer Liste, nach ihrer Giite sortiert und die n% schlechtesten bzw. die n schlechtesten MeBdaten aussortiert wer- 
den. Insbesondere ist zu priifen, ob die n% bzw. die n schlechtesten MeBdaten in einem zusammenhangenden Bereich 

60 liegen. Ist dies der Fall, so werden diese MeBdaten nicht aussortiert, da sie mit hoher Wahrscheinlichkeit keine MeBfeh- 
ler, sondern einen zusammenhangenden Bereich bestimmen, der durch das Ersatzmodell nicht ausreichend genau abge- 
bildet wurde. 

Eine andere Weiterbildung besteht darin, daB die MeBdaten einer Vorverarbeitung unterzogen werden. Da in einem 
realen vorgegebenes System eine groBe Menge an MeBdaten pro Zeiteinheit anfallen, ist es sinnvoll, diese MeBdaten ei- 
65 ner Vorverarbeitung zu unterziehen und somit zu gewahrleisten, daB weitgehend signifikante MeBdaten in die Bildung 
des Ersatzmodells einflieBen. Bevorzugt findet die Vorverarbeitung in einer Reduktion der Anzahl von MeBdaten ihre 
Auspragung. 

Dabei werden die MeBdaten entsprechend vorgegebener Kriterien in Klassen eingeteilt. Die MeBwerte einer Klasse 
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en MeBwerte, deren Bewertung unterhalb eines vorgegebSme 



werden bewertet UDd dieJH^en MeBwerte, deren Bewertung unterhalb eines vorgegebelwersten Schwellwertes liegt, 
werden aussortiert. Durch die Aussortierung der MeBwerte ergibt sich eine Reduktion hinsichdich der Anzahl der MeB- 
werte. Somit liegen fur eine Weiterverarbeitung eine deutlich reduzierte Anzahl von MeBwerten vor. Die Weiterverarbei- 
tung kann mit gegeniiber der nicht reduzierten Anzahl von MeBwerten geringerem Rechenaufwand erfolgen. 

Auch konnen die Klassen selbst bewertet werden. Insbesondere kann eine Klasse, deren Bewertung unterhalb eines 5 
vorgegebenen zweiten Schwellwerts liegt, aussortiert werden. Hierdurch ergibt sich eine zusatzliche Reduktion der An- 
zahl der MeBwerte. 

Eine andere Weiterbildung der Vorverarbeitung besteht darin, daB ein Kriterium fur die Klasseneinteilung darin be- 
steht, daB pro Klasse MeBwerte zu einer Vorgabe von Einstellparametern des technischen Systems bestimmt werden. Ty- 
pischerweise wird das technische System anhand einer vorgegebenen Anzahl von Einstellparametern eingestellt, nach 10 
Einstellung erfolgt eine (zumeist zeitverzogerte) Reaktion des Systems auf die Einstellparameter (Einschwingverhalten, 
Einschwingvorgang des Systems) Nach Einstellung werden somit eine bestimmte Menge an MeBwerten aufgenommen, 
die dem Einschwingvorgang zugeordnet werden konnen, wobei nach abgeschlossenem Einschwingvorgang (Ubergang 
in den stationaren Betrieb) weiterhin MeBwerte anfallen, die dem vorgegebenen Satz Einstellparameter zugeordnet wer- 
den. Durch Verstellung der Einstellparameter wird eine neue Klasse bestimmt. Alle MeBwerte, die jeweils nach einer 15 
Verstellung der Einstellparameter anfallen, gehoren in eine eigene Klasse. 

Zusatzlich konnen MeBwerte einer Klasse, die dem jeweiligen Einschwingvorgang zuordenbar sind, aussortiert wer- 
den. 

Weiterhin konnen fehlerhafte MeBwerte aussortiert werden. Die Einstellung groBer technischer Systeme ist in vielen 
Fallen auf eine langfristigen stationaren Betrieb ausgerichtet. MeBwerte, die sich auf den Einschwingvorgang (von kur- 20 
zer Dauer im Verhaltnis zum stationaren Betrieb nach abgeschlossenem Einschwingvorgang) beziehen, werden sinnvoll 
aussortiert, da durch sie MeBwerte fur den stationaren Betrieb verfalscht werden. Insbesondere im Rahmen einer Model- 
lierung des technischen Systems, sind die MeBdaten des stationaren Verhaltens des technischen Systems interessant. 

Eine Ausgestaltung besteht darin, die Anzahl der MeBwerte in einer Klasse dadurch zu reduzieren, daB mindestens ein 
reprasentativer Wert fur die MeBwerte der Klasse bestimmt wird. Solch ein reprasentativer Wert kann sein: 25 

a) ein Mittelwert (z. B. ein gleitender Mittelwert) der MeBwerte der Klasse, 

b) ein Maximalwert der MeBwerte der Klasse, 

c) ein Minimal wert der MeBwerte der Klasse, 

d) ein Median. 30 

Bei Variante d) liegt ein Vorteil darin, daB immer ein Wert bestimmbar ist, den es tatsachlich gibt, wohingegen der Mit- 
telwert a) selbst nicht als Wert vorkommt. 

Je nach Anwendungsfall, kann eine geeignete Wahl zur Bestimmung des reprasentativen Werts einer Klasse erfolgen. 

Eine ganze Klasse mit MeBwerten kann aussortiert werden, wenn diese weniger als eine vorgegebene Anzahl MeB- 35 
werte enthalt. 

Eine andere Ausgestaltung besteht darin, daB diejenigen MeBwerte aussortiert werden, die um mehr als einen vorge- 
gebenen Schwellwert von einem vorgebbaren Wert verschieden sind. Der vorgebbare Wert kann ein Mittelwert aller 
MeBwerte der Klasse oder ein zu erwartender MeBwert auf die jeweiligen Einstellparameter des technischen Systems 
sein. 40 

Im Rahmen einer anderen Weiterbildung werden die mittels Entwurf gewonnenen Daten zur Steuerung einer techni- 
schen Anlage eingesetzt. Zusatzlich kann die Steuerung der technischen Anlage zur Laufzeit des Systems, also Online, 
erfolgen. 

Auch wird zur Losung der Aufgabe eine Anordnung zum Entwurf eines technischen Systems angegeben, die eine Pro- 
zessoreinheit aufweist, welche Prozessoreinheit derart eingerichtet ist, daB MeBdaten eines vorgegebenen Systems an- 45 
hand eines Ersatzmodells beschreibbar sind. Eine MaBzahl fur die Giite des Ersatzmodells ist ermittelbar, indem die 
MeBdaten des vorgegebenen Systems mit durch das Ersatzmodell bestimmten Daten verglichen werden. Aus der MaB- 
zahl fur die Giite ist das Ersatzmodell dahingehend anpaBbar, daB es eine moglichst hohe Giite aufweist. Das hinsichtlich 
seiner Giite angepaBte Ersatzmodell ist zum Entwurf des technischen Systems einsetzbar. 

Diese Anordnung ist insbesondere geeignet zur Durchfuhrung des erfindungsgemaBen Verfahrens oder einer seiner 50 
vorstehend erlauterten Weiterbildungen. 

Ausfuhrungsbeispiele der Erfindung werden nachfolgend anhand der Zeichnung dargestellt und erlautert. 

Es zeigen 

Fig. 1 ein Blockdiagramm, das Schritte eines Verfahrens zum Entwurf eines technischen Systems enthalt; 

Fig. 2 eine schematische Skizze eines Recovery-Boilers; 55 

Fig. 3-5 EingangsgroBen, StellgroBen und AusgangsgroBen des Recovery-Boilers. 

In Fig. 1 ist ein Blockdiagramm dargestellt, das Schritte eines Verfahrens zum Entwurf eines technischen Systems ent- 
halt. In einem Schritt 101 wird anhand von MeBdaten ein Ersatzmodell gebildet. Bevorzugt ist dieses Ersatzmodell ein 
Regressionsmodell. Um das in Schritt 101 entstandene Ersatzmodell auf die MeBdaten anzupassen, also eine Verfeine- 
rung des Ersatzmodells vorzunehmen, so daB die MeBdaten das Ersatzmodell in ausreichender Naherung beschreiben, 60 
wird in einem Schritt 102 eine MaBzahl fur die Giite des Ersatzmodells bestimmt. Diese MaBzahl wird bestimmt, indem 
die MeBdaten des vorgegebenen Systems mit durch das Ersatzmodell bestimmten Daten verglichen werden. Bevorzugt 
erhalt jedes MeBdatum eine MaBzahl fur die Giite, welche MaBzahl die Abweichung des MeBdatums mit dem zugehori- 
gen von dem Ersatzmodell bestimmten Wert kennzeichnet Die Summe aller MaBzahlen fur die Giite fur alle MeBdaten 
bestimmt eine Gesamtgiite fur das Ersatzmodell. In einem Schritt 103 wird die Giite maximiert, indem die MaBzahl fur 65 
die Giite (bzw. eine negative Giite fur die Ubereinstimmung des Ersatzmodells mit dem vorgegebenen System) minimiert 
wird. Ist eine entsprechend hohe Giite fur das Ersatzmodell bestimmt, so wird in einem Schritt 104 dieses Ersatzmodell 
fur den Entwurf des technischen Systems eingesetzt Der Entwurf kann sowohl ein Neuentwurf (vgl. Schritt 105) oder 
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eine Anpassung eines schc^^^ehenden technischen Systems (vgl. Schritt 106) sein. 

Fig. 2 zeigt eine schematische Skizze eines Recovery-Boilers. Nachfolgend wird anhanddes Bei spiels "Recovery- 
Boiler" ein Ausfuhrungsbeispiel des oben beschriebenen Verfahrens veranschaulicht. 

In der Papier- und Zellstoffindustrie werden zum AufschluB von Zellstoff verschiedene Chemikalien sowie Warme 
5 und Elektroenergie bendtigt. Aus einer eingedeckten ProzeBablauge (Schwarzlauge) lassen sich mit Hilfe dcs Recovery- 
Boilers die verwendeten Chemikalien und zusatzlich Warmeenergie zuriickgewinnen. Ein Grad fur die Zuriickgewin- 
nung der Chemikalien ist von entscheidender Bedeutung fur die Wirtschaftlichkeit der Gesamtanlage. 

Die Schwarzlauge wird in einem Schmelzbett 201 verbrannt. Dabei bildet sich eine Alkalischmelze, die iiber eine Lei- 
tung 202 abflieBt. Aus den Bestandteilen der Alkalischmelze werden in weiteren Verfahrensschritten die eingesetzten 
10 Chemikalien zuriickgewonnen. Frei gewordene Verbrennungs warme wird zur Erzeugung von Wasserdampf genutzt. Die 
Verbrennung der Ablauge und damit die Ruckgewinnung der Chemikalien beginnt mit der Zerstaubung der Schwarz- 
lauge iiber Zerstauberdusen 204 in eine Brennkammer 203. Partikel der zerstaubten Schwarzlauge werden bei ihrem Fall 
durch das heiBe Rauchgas getrocknet. Die getrockneten Laugenpartikel fallen auf das Schmelzbett 201, wobei eine erste 
Verbrennung und eine chemische Reduktion stattfinden. Fluchtige Bestandteile und Reaktionsprodukte gelangen in eine 
15 Oxidationszone, in der oxidierende Reaktionen ablaufen und in der die Verbrennung abgeschlossen wird. 

Wichtige Zielvorgaben fiir die Steuerung des Recovery-Boilers sind die Dampfproduktion zur Energiegewinnung, die 
Einhaltung von Emissions werten unter Umweltgesichlspunkten und die Effizienz der chemischen Reduktion. 

Der Verbrennungsvorgang, und damit die Zielvorgaben, werden insbesondere durch die Luftzufuhr in drei Ebenen 
(Primary Air (PA), Secondary Air (SA), Tertiary Air (TA)) gesteuert. Der GesamtprozeB unterliegt zahlreichen Einfliis- 
20 sen, die bei der Modellierung zu beriicksichtigcn sind: 

a) Die Messung der GroBen unterliegen oftmals starken Schwankungen; 

b) Es existieren nicht gemessene und nicht meBbare EinfluBgroBen; 

c) Bei jeder Anderung der einstellbaren Parameter kornmt es zu Einschwingvorgangen; 

25 d) Die technische Anlage verschmutzt und wird in vorgegebenen Abstanden gereinigt, wodurch im Systemverhal- 

ten jeweils ein zeitlicher Drift bewirkt wird. 

Die gemessenen GroBen des Gesamtprozesses werden in EingangsgroBen (vgl. Fig. 3) und AusgangsgroBen (vgl. Fig. 
5) unterteiit. Jede Minute werden MeBwerte abgespeichert. Vier der EingangsgroBen sind gleichzeitig auch StellgroBen 

30 (auch: einstellbare Parameter; vgl. Fig. 4). Die StellgroBen sind im wesentlichen als unabhangig voneinander einstellbare 
freie Parameter des Gesamtprozesses anzusehen. Einige der anderen EingangsgroBen sind von den StellgroBen mehr 
oder minder abhangig. GemaB einer Vorgabe sind beim Recovery-Boiler die GroBen "BL Front Pressure" und "BL Back 
Pressure" stets gleich zu regeln. Die vier StellgroBen (vgl. Fig. 4) sind vorzugsweise abzuspeichern ais StellgroBen (mit 
dem gewiinschten, voreingestellten Wert) und als EingangsgroBen (mit dem gemessenen, realen Wert). 

35 Beim Recovery-Boiler besteht eine Problemstellung darin, in Abhangigkeit von den einstellbaren Parametem be- 
stimmte Zielvorgaben, die iiber gemessene GroBen definiert werden, zu erfullen. Hier wird eine dreistufige \fcrgehens- 
weise zur Losung des Problems gewahlt: 

1. Die zu betrachtenden Zielvorgaben werden durch stochastische Methoden modelliert, wobei diese Modelle 
40 durch neue Messungen aktualisiert werden (datengetriebene, empirische Modellierung). Dabei ist es sinnvoll, nicht 

nur ein einziges Modell zu verwenden, sondem globale Modelle fiir die Identifikation interessanter Gebiete in ei- 
nem durch die Zielvorgaben bestimmten Parameterraum und lokale Modelle zur exakten Berechnung optimaler Ar- 
beitspunkte einzusetzen. Die verwendeten Modelle werden durch Gutemafle bewertet. 

2. Falls die betrachteten Modelle aufgrund der Datenlage nicht hinreichend genau sind (GiitemaBe), werden gezielt 
45 neue Arbeitspunkte zur Modellverbesserung ausgewertet (Experimental Design). Ferner werden durch Verwendung 

globaler stochastischer Optimierungsverfahren bzgl. der Zielvorgaben attraktive Gebiete in Abhangigkeit vom ak- 
tuellen globalen Modell identiflziert. 

3. Fiir die lokale Optimierung werden lokale Modelle konstruiert und die zur Verfugung stehenden Datensatze ge- 
gebenenfalls gezielt erweitert (Experimental Design). 

50 

Bei den Zielvorgaben handelt es sich um physikalischtechnische bzw. betriebswirtschaftliche Kriterien, die in der Re- 
gel Randbedingungen und/oder Sicherheitsbedingungen entsprechen mussen. Haufig sind mehrere dieser Kriterien 
gleichzeitig zu betrachten. Die Verwendung eines stochastischen Modells kann insbesondere dazu verwendet werden, 
die zu optimierenden ZielgroBen und ihre Abhangigkeit von den einzustellenden Parametern im Rechner zu simulieren. 

55 Dies ist dann notwendig, wenn Messungen sehr kostenintensiv bzw. sehr zeitaufwendig sind. Bei Sicherheitsanforderun- 
gen konnen mcgiichc Gefahrensituationen vennieden werden. 

Beim Recovery-Boiler ist eine Online-Optimierung, die auf mehreren Daten basiert, notwendig, weil die physikalisch- 
chemischen Prozesse nicht mit ausreichender Genauigkeit quantitativ modelliert werden konnen und weil das ^Vbrhalten 
der Anlage im Verlauf des Betriebs Schwankungen unterliegt. Das Wissen iiber dieses Verhalten muB stetig durch ge- 

60 zielte Wahl neuer Arbeitspunkte erweitert werden. Daher empfiehlt sich im Rahmen der Online-Optimierung das bereits 
beschriebene dreistufige Vorgehen der stochastischen Modellierung und der mathematischen Optimierung. 

BESCHREffiUNG DER EINGANGSGROSSEN 

65 Die a EingangsgroBen (a e N, N: Menge der natiirlichen Zahlen) sind im allgemeinen von n StellgroBen n e N und 
von Zufallseffekten abhangig. Sie konnen wie folgt beschrieben werden: 

Es seien (ft S, P) ein Wahrscheinlichkeitsraum und B v eine Borelsche 0-Algebra iiber R v (R: Menge der reellen Zahlen) 
fiir jedes v e N. Die EingangsgroBen werden iiber eine B n X S - B a -meBbare Abbildung (p dargestellt: 
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<p:R n xQ — R a (1).^^ 

Die Definitionsmenge der Abbildung j ist ein kartesisches Produkt zweier Mengen. Betrachtet man die jeweiligen Pro- 
jektionen auf die Einzelmengen, so erhalt man folgende Abbildungen: 

5 

<p x : Q. — R a ,C0 — <p(x,CD) fur alle x e R n (2), 
<p to :R D -.R a ,x-^(p(x,(0)furaile(a) e n (3). 

{(px; x e R n } ist ein stochastischer ProzeB mit einer Indexmenge R D und eine Abbildung (j/ 0 ist fur jedes Ereignis (06 10 
Q ein Pfad dieses stochastischen Prozesses. 

Beim Recovery-Boiler ist n = 4 und a = 14 (nach Elimination der GroBe "BL Back Pressure"). 

Aufgrund der geforderten MeBbarkeit der Abbildung (f> x ist fur jedes x e R n die Abbildung <p x eine Zufallsvariable. 
Unter geeigneten Zusatzvoraussetzungen konnen Erwartungswerte und hohere Momente betrachtet werden. Dieser Zu- 
gang ermoglicht den Schritt von stochastischen Modellen zu deterministischen Optimierungsproblemen. Bei einem de- 15 
terministischen Optimierungsproblem ist die Zielfunktion direkt mittels einer Variablen einstellbar, wohingegen die sto- 
chastische GroBe die Zielfunktion beeinfluBt, aber keine gezielte Einstellung ermoglicht. 
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BESCHRETOUNG DER AUSGANGSGROSSEN 

Das ProzeBmodell M des Recovery-Boilers wird als Funktion in Abhangigkeit von den EingangsgroBen und weiteren 
ZufallsefFekten beschrieben. Dabei sei (Q S, P) der obige Wahrscheinlichkeitsraum. Das ProzeBmodell M ist dann eine 
B a x S - B b -meBbare Abbildung: 

M:R a x£2^R b (4), 25 

wobei b die Anzahl der AusgangsgroBen bezeichnet. 

Da der Recovery-Boiler einem zyklischen zeitlichen Drift unterliegt (von Reinigungsphase zu Reinigungsphase), ist 
zudem eine Beschreibung mit einem Zeitparameter denkbar. Die AusgangsgroBen lassen sich durch B n x S - B -meBbare 
Abbildungen \j/ darstellen: 30 

\|/:R n xn — R b (5), 

(x,o>)^M((p(/x,G)),G)) (6). 

35 

Betrachtet man die jeweiligen Projektionen auf die Einzelmengen der Definitionsmenge, so erhalt man folgende Ab- 
bildungen 

\j/ x : Q— ► R b ,w -+ \|/(x,co) fur alle x e R n (7), 

40 

y 05 : R n — * R b ,x — * \y(x,co) fur alle (0€Q (8). 

{Xj^; x e R n } ist einstochastischer ProzeB mit einer Indexmenge R n und die Abbildung \|/*° ist fur jedes fi> e Q ein 

Pfad dieses stochastischen Prozesses. 

Beim Recovery-Boiler ist b = 15. 45 
Die Tatsache, daB bei der Definition von y zwischen den verwendeten Ereignissen CO nicht unterschieden wird, bedeu- 

tet keine Einschrankung, da Q. als kartesisches Produkt aus einem Qi und einem Q.2 dargestellt werden kann. Die obige 

Darstellung umfaBt somit auch das Modell: 

\|f:R'xfl 1 xQ 2 ^R b (9), 50 



(x,©!,^) -> MCq^x,©!),^) (10). 

BESCHRETOUNG DER ZUR VERFUGUNG STEHENDEN DATENSATZE 

55 

Mit den Beschreibungen in den beiden vorangegangenen Abschnitten kann man die EingangsgroBen und die Aus- 
gangsgroBen gemeinsam zu MeBgroBen (= MeBdaten) <P zusammenfassen. Oist eine B n x S - B m -meBbare Abbildung 
mit m = a + b und 

0 : R n x £2 — ► R m (11), 60 




(12) . 
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Betrachtet man wieder die jeweiligen Projektionen auf die Einzelmengen der Definitionsmenge, so erhalt man fol- 
gende Abbildungen: 
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<D X : Q — R ra ,G) — O(x,(0)^PBe x e R n (13), 

O w : R n — R m ,x — <I>(x,G)) fur aUe <o e Q (14). 

5 {O x ; x e R D ) ist ein stochastischer ProzeB mit einer Indexmenge R n und die Abbildung <J> W ist fiir jedes G> e ft ein 
Pfad dieses stochastischen Prozesses. 

Fiir jedes gewahlte StellgroBentupel x werden beim Recovery-Boiler viele Realisiemngen von <& x ermittelt und abge- 
speichert, d. h. zu jedem Xj e R n werden zahlreiche Realisiemngen 

10 <Dj k :=*Kxj,(Dj k ) (15) 
mit 

G>j k eii;k=l, 2,. . .,vj; 
15 Vj e N; j= 1, 2,. . .,u; u e N 

betrachtet. Die gespeicherten Datensatze Dj k des Recovery-Boilers sind also (n + m)-Tupel: 



20 D i k = I /i>j k J ' k = • • • r v j ; j = 1,2, . . . f u (16) 

Dabei wird Dj lkl vor Dj 2k2 abgespeichert, wenn 
25 Qi<j2) V (0i=j2) A (k!<k 2 )) 
gilt. 

DATENKOMPRES SION DURCH KLAS SENEINTEELUNG DER PARAMETER 

30 

Da fiir jedes StellgroBentupel x im allgemeinen mehrere Realisiemngen von O x vorliegen, bietet sich aufgrund der 
komplexen stochastischen Eigenschaften des zu betrachtenden Prozesses als erster Schritt der statistischen Datenanalyse 
eine Klasseneinteilung der Parameter durch Bildung arithmetischer Mittelwerte an. Zudem werden offensichtlich fehler- 
hafte Datensatze ausgesondert. Ein offensichtlich fehlerhafter Datensatz ist bspw. eine physikalisch unmogliche Mes- 
35 sung, der insbesondere aufgmnd einer vorgenommenen Einstellung real gar nicht vorkommen kann. 

Vorgehensweise: 

1. Datensatze, bei denen die GroBe "BL Front Pressure" ungleich der GroBe "BL Back Pressure" ist, werden aus- 
sortiert, da diese beiden Werte nach Vorgabe der Anlagensteuerung gleich sein mussen. Der Datenverlust ist sehr 

40 gcring. 

2. Die Datensatze werden auf Klassen aufgeteilt, in denen die vier Einstellparameter (PA, S A, TA, BL Front Pres- 
sure, siehe oben) zeitlich aufeinanderfolgend konstant sind, d. h. die j-te Klasse besteht aus den Datensatzen Dj ♦ . 

3. Klassen, in denen sich weniger als 30 Datensatze befinden, werden ausgesondert, damit Einschwingvorgange 
keinen groBen EinfluB haben. 

45 4. Fiir jede Klasse werden ein arithmetischer Mittelwert <J>j und eine empirische Standardabweichung Sj fiir alle 

MeBgroBen ermittelt: 

l" Vj 

°3 = —■ I>jk (17), 

so V : k=i 
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S-i = 



Vj-1 



Vj-1 



I Hi 

k=l 



k=l 



v 



(18) . 



5. Klassen bei denen die Mittelwerte fiir die GroBen PA, SA, TA oder BL Front Pressure zu weit von den entsprechen- 
den Einstellparametem entfernt sind, werden ausgesondert. Auf diesen Klassen konnten also die Einstellwerte nicht er- 
reicht werden. 
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IEN FUR DIE GEGEBENEN KLASSEN UND IHRE G: 



CHE DARSTELLUNG 



Neben den arithmetischen Mittelwerten und den empirischen Standardabweichungen, die fiir die einzelnen Klassen 
bestimmt wurden, wird noch eine gemeinsame Standardabweichung s bestimmt gemaB 



S = 



V - 1 



v - 1 



- *f 



- l 



^2 



(19) 



10 



15 



Dabei steht u fiir die Anzahl der Klassen (hier 205) und v fur die Summe der Vj, d. h. v ist die Anzahl aller verwendeten 
MeBwerte (hier 389 15). 20 

LINE ARE REGRESSIONSMODELLE FUR FUNKTIONSAPPROXIMAnONEN 

Fiir jede MeBgroBe (MeBdatum) (i = 1, 2,. . .,m) wird, basierend auf dem arithmetischen Mittel iiber die Klassen, 
ein lineares Regressionsmodell in Abhangigkeit von der quadratischen Kombination der vier Einstellparameter berech- 25 
net. In der folgenden Darstellung ist x E R 4 , wobei 
x (l) : Primary Air (PA) 
x (2) : Secondary Air (SA) 
x (3) : Tertiary Air (TA) 

x (4) : Black Liquor (BL) Front Pressure 30 
gilt, u e N bezeichnet die Anzahl der Klassen. Jede MeBgroBe wird durch 

0^(x,<d) = a[r(x) + ei(co) (20) 

mit ai e R 15 modelliert. Dabei gelten 
r:R 4 — R 15 (21) 

feu fe' ^3/ ^4) T "> ( 1# Clr £>2, £>3> l-A* & & & & 

Cl&rM3.ClC4.C2C3.«4.C3?4) T (22), 
d. h. Polynome zweiten Grades werden an die MeBdaten angepaBt, und 

e^O^R (23) 

ist eine Zufallsvariable mit Erwartungswert 0. 

Der Vektor a 4 wird mit der Methode der kleinsten Fehlerquadrate bestimmt, allerdings werden statt der Originaldaten- 
satze 



35 



40 



45 



die arithmetischen Mittel 

verwendet. Diese Vorgehensweise eignet sich, da durch lineare Regressionsmodelle insbesondere Erwartungswerte ge- 
schatzt werden. Somit erhalt man das Minimierungsproblem: 



min 



15 



0<i) 



r(xi) T 
< r ( x uF; 



(15) 



2 J 



(24) 
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Es sei aj der optimale Pil^Pes quadratischen Minimierungsproblems aus Gleicbung (S^Rrner gelte: 



yi : = 









\*i J 



G R 



u 



(25), 



10 



15 



25 



U 



7 . := ±. Z o(i) eR 
j=i 



(26) 



Zur Validierung des Regressionsansatzes wird ein BestimmtheitsmaB R 2 berechnet gemaB 

v2 



20 R 2 



mit 



rp O 

2 ^(i)T^(i) _ „^7? 



^(i) _ W - uy 



(27) 



30 



35 



40 



45 



50 



55 



60 



65 



*0) = 



(28) . 



Je naher Rf bei 1 liegt, desto besser wird die abhangige Variable durch die unabhangigen Variablen dargestellt 
(0 < Rf < 1). 

Zusatzlich wird ein Maximum EJjU fur einen Absolutwertes der Abweichung der Daten vom Model! angegeben als 



iax : = max \ 
j = l,...,u I 



(29) 



Ego% derjenige Wert, unterhalb dessen mindestens 90% der Absolutwerte der Abweichungen der Daten vom Modell 
liegen. Analog dazu ist derjenige Wert, unterhalbdessen mindestens 80% der Absolutwerte der Abweichungen der 
Daten vom Modell liegen. Mit dem optimalen Punkt a[ des Minimierungsproblems gemaB Gleichung (24) laBt sich ein 
Modell<l> (l) des Erwartungswertes der MeBgroBe angeben zu 

0> a) : = R n — R (30), 

x — * T r(x) (31). 

Insbesondere laBt sich der Gradient VU> (l) analytisch angeben mit 

fur alle x e R n 



= — (x) 
dx 



a-i 



(32) 



Patentanspriiche 



1. Verfahren zum Entwurf eines technischen Systems, 

a) bei dem MeBdaten eines vorgegebenen Systems anhand eines Ersatzmodells beschrieben werden; 

b) bei dem eine MaBzahl fur die Gute des Ersatzmodells ermittelt wird, indem die MeBdaten des vorgegebe- 
nen Systems mit durch das Ersatzmodell bestimmten Daten verg lichen werden; 

c) bei dem aus der MaBzahl fur die Gute das Ersatzmodell dahingehend angepaBt wird, daB es eine moglichst 
hohe Gute aufweist; 
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d) bei dern dSBKichtlich seiner Giite angepaBte Ersatzmodell zum Entwurr^Kechmscheri Systems einge- 
setzt wird. 

2. Verfahren nach Anspruch 1, bei dem das Ersatzmodell ein Regressionsmodell ist. 

3. Verfahren nach Anspruch 1 oder 2, bei dem die Giite anhand einer quadratischen Abweichung der MeBdaten von 
den durch das Ersatzmodell bestimmten Daten ermittelt wird. 5 

4. Verfahren nach einem der vorhergehenden Anspruche, bei dem die MeBdaten nach deren Giite, bezogen auf de- 
ren Abweichung von den durch das Ersatzmodell bestimmten Daten, sortiert werden und eine vorgegeben Anzahl 
von n% schlechtesten MeBdaten aussortiert werden. 

5. Verfahren nach Anspruch 4, bei dem die n% schlechtesten MeBdaten nicht aussortiert werden, wenn sie in einem 
zusammenhangenden Bereich liegen. 10 

6. Verfahren nach einem der vorhergehenden Anspruche, bei dem die Menge der MeBdaten im Rahmen einer \fc>r- 
verarbeitung reduziert wird. 

7. Verfahren nach Anspruch 6, bei dem die Vorverarbeitung eine Klasseneinteilung der MeBdaten umfaBt. 

8. Verfahren nach einem der vorhergehenden Anspruche, bei dem die mittels Entwurf gewonnenen Daten zur 
Steuerung einer technischen Anlage eingesetzt werden. 15 

9. Verfahren nach Anspruch 8, zur Online- Anpassung der Steuerung fur die technische Anlage. 

10. Anordnung zum Entwurf eines technischen Systems, mit einer Prozessoreinheit, die derart eingerichtet ist, daB 

a) MeBdaten eines vorgegebenen Systems anhand eines Ersatzmodells beschreibbar sind; 

b) eine MaBzahl fur die Giite des Ersatzmodells ermittelbar ist, indem die MeBdaten des vorgegebenen Sy- 
stems mit durch das Ersatzmodell bestimmten Daten verglichen werden; 20 

c) aus der MaBzahl fur die Giite das Ersatzmodell dahingehend anpaBbar ist, daB es eine moglichst hohe Giite 
aufweist; 

d) das hinsichtlich seiner Giite angepaBte Ersatzmodell zum Entwurf des technischen Systems einsetzbar ist. 



Hierzu 3 Seite(n) Zeichnungen 25 
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